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ABSTRACT 

We are engaged in a comprehensive program to find reliable elemental abun- 
dances in and to probe the physical structure of the Orion Nebula, the brightest 
and best-resolved H II region. In the course of developing a robust extinction 
correction covering our optical and ultraviolet FOS and STIS observations, we 
examined the decrement within various series of He I lines. The decrements of 
the 2 3 S — n 3 P, 2 3 P — n 3 S and 3 3 S — n 3 P series are not in accord with 
case B recombination theory. None of these anomalous He I decrements can be 
explained by extinction, indicating the presence of additional radiative transfer 
effects in He I lines ranging from the near-IR to the near-UV. CLOUDY photoion- 
ization equilibrium models including radiative transfer are developed to predict 
the observed He I decrements and the quantitative agreement is quite remark- 
able. Following from these results, select He I lines are combined with H I and 
[O n] lines and stellar extinction data to validate a new normalizable analytic 
expression for the wavelength dependence of the extinction. In so doing, the 
He + /H + abundance is also derived. 

Subject headings: H II regions — ISM: individual (Orion Nebula) — ISM: abun- 
dances — dust, extinction 
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1. Introduction 

The Orion Nebula is the brightest and best-resolved H II region and is the defining case 
with a blister geometry (Zuckerman 1973; Balick et al. 1974). Our study of this region makes 
use of emission-line diagnostics from HST FOS and STIS spectra supplemented by extensive 
ground-based observations (e.g., Baldwin et al. 1996; Rubin et al. 1997; Baldwin et al. 2000), 
with the goal of probing the nebula's physical structure and determining a robust set of its 
elemental abundances. 

To complete a full spectral analysis in the UV and visible, we have had to develop 
and test an extinction curve amended slightly from those often applied to visible Orion 
Nebula observations (Costero & Peimbert 1970; Cardelli et al. 1989, hereafter CP70 and 
CCM89, respectively). See § 3. In the course of this work we examined the decrement within 
various series of He I lines, planning to use these lines to constrain extinction, primarily in 
the ultraviolet. As reported in Martin et al. (1996) and now fully described in § 4, our 
FOS and STIS observations have revealed a decrement within the ultraviolet 2 3 S — n 3 P 
series that is not in accord with case B (as defined in Baker & Menzel 1938) recombination 
theory (Smits 1996; Benjamin et al. 1999; Porter et al. 2005), and cannot be explained by 
extinction effects. This anomalous decrement stems from the metastability of 2 3 S, leading 
to radiative transfer effects which then affect other series (Robbins 1968; Osterbrock 1989). 
We compare the predictions of this theory with our set of HST FOS and STIS observations 
- the first detailed observational analysis of He I UV lines originating from high n terms. In 
a comprehensive examination including other datasets we also show quantitatively how the 
same theory self-consistently accounts for resonance fluorescence enhancement observed in 
several other lines of two related series. Including the radiative transfer effects reduces \ 2 
by a factor of roughly 10 as compared to case B recombination alone, showing a remarkable 
agreement between radiative transfer theory and observations. 



2. Observations 

HST FOS and STIS spectra were obtained with spectral coverage from the UV (1600A) 
through to the visible (7400A) for two lines of sight (roughly 1SW and x2, Baldwin et al. 
1996; Rubin et al. 1997). The FOS aperture has a diameter 0'.'86 and the STIS slit is 52" x 0'.'5 
(refer to Fig. 1) of which only (the roughly central) 28" are recorded using the STIS UV 
detectors (FUV-MAMA, NUV-MAMA) . The UV and optical spectra were adjusted to align 
spatially using a bright feature common to all spectra: proplyd 159-350 (O'Dell & Wen 
1994). 
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The placement of Slit 1 (see Fig. 1) was such that it was centered on a region of the 
nebula where there is an extinction gradient. Because of this, we divided the 52" x 0"5 slit 
into four 13" x &!5 sections 1 , labeled "a" through "d". The different spatial coverage of the 
MAMA and CCD detectors means that the outer two sections (SLITla and SLITld) have 
reliable coverage only in the visible, but the central two sections (SLITlb and SLITlc) have 
coverage in the visible and the UV. The same sectioning was done for Slit 2 (x2, see Fig. 1). 
However, there was little change in the extinction across Slit 2 so the two central sections 
(SLIT2b and SLIT2c) should be similar. 

The reduction of STIS spectra was performed as in Rubin et al. (2003) and the reduction 
of FOS spectra as in Rubin et al. (1998). Our HST STIS (SLITlb, SLITlc, SLIT2b, SLIT2c) 
and FOS (1SW) observations are presented in Table 1. 

A wavelength-dependent extinction correction must still be applied to the entire spec- 
trum (UV through visible/near-IR). To date, there has not been a robust way of reliably 
correcting such a large portion of the spectrum to allow for full (i.e., UV through visible) 
spectral analysis. However, with the broad spectral coverage of the HST data, we are able to 
develop and test a smooth consistent extinction correction curve. The development and vali- 
dation of such an analytic extinction curve is discussed in § 3. To ensure that any extinction 
curve is well constrained, we must maximize the number of lines used in its determination. 
In § 4, we discuss the basis for the inclusion of He I (especially He I UV) lines. 



3. Extinction corrections 

Were it not for extinction, the Orion Nebula would appear much brighter. Typical 
optical depths reported are 7y ~ 1 . From the shape of the extinction curve (e.g., see curves 
in Kim et al. 1994), the optical depth in the Lyman continuum (T912) would be several 
times this value. The relatively soft radiation from the primary ionization source (O star, 
9 1 Ori C) cannot photoionize an optical depth to dust of much more than rg 12 ~ 1 indicating 
that a large part of the optical extinction must come from neutral material in the foreground. 
Further evidence for and properties of this "veil" are summarized by Abel et al. (2004). 

All observed emission-line surface brightnesses or fluxes, from the infrared to ultra- 
violet, are affected by this wavelength-dependent extinction and must be corrected before 



1 All STIS spectra have small fiducial bar sections which interfere with the spatial coverage of the spectra. 
For this reason, when average surface brightnesses of tiles are being calculated along the slit, these fiducial 
bar sections are avoided. 
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further analysis and interpretation. Extinction is characterized by its wavelength dependence 
(shape) and the amount at a given wavelength (amplitude). 

Stellar extinction deduced from stars in the nebula provides a useful reconnaissance of 
the shape of the Orion Nebula extinction curve, providing a continuous curve with which 
to make predictions (by interpolation) for any given line observed. We have reconsidered 
the wavelength dependence of the infrared, optical and ultraviolet extinction, presenting 
our results as a convenient analytic expression describing this shape ( "normalized extinction 
curve") in § 3.1.1. The amplitude by which to scale the normalized extinction curve can be 
determined for a given line of sight from nebular lines whose ratios of line emissivities are 
known from atomic theory (e.g., H I Balmer series). An optimal implementation is discussed 
in § 3.1.2. In § 3.2, we compare our extinction curve with the nebular extinction curve of 
CP70. In § 3.3, we discuss our extensive validation of the use of our stellar-based extinction 
curve for interpolating corrections to the nebular emission in the optical, the ultraviolet and 
the near-infrared. 



3.1. Parameterization of extinction 

We have already stressed how one requires a shape (s\) and amplitude (a) to describe 
the extinction: 

T X = as x . (1) 

For different lines of sight, grain properties like size distribution and composition are mani- 
fested in s\, while a scales with the column density of dust. 

Let F represent the observed flux of a patch of nebula (surface brightness times solid 
angle) or of a star, with subscript o indicating the intrinsic flux - as it would have been had it 
not been affected by extinction r. (The dependence on wavelength A has been suppressed.) 
Thus F/F Q = exp(— r). For a uniformly bright nebular patch this is also the ratio of the 
specific intensities, I jl . Note that I can be predicted by atomic theory, from the emissivity 
j and (for most lines) the emission measure (Osterbrock 1989). We adopt I Q oc j in the 
applications here, so that 

Fa = constant x j\ exp(— as\). (2) 

The logarithmic (base 10) extinction is C = — \og(F / F Q ) = — log(exp(— r)) = 0.4343t. 
The extinction in magnitudes is A = — 2.5 log(F/F ) = 2.5C = 1.086r. Given the linear 
transformations between A, C, and r, a ratio like A\ 1 /A\ 2 evaluated for two wavelengths is 
clearly interchangeable with the same ratio for C or r. 
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In the literature there are various parameterizations of A, C, and r that are worth 
distinguishing, since they imply different choices and normalizations of s\. Osterbrock (1989) 
uses S\ = 0.5t a /(th 7 — r Ha ) and comments on the utility of the differential form s x — s H( 3- 
Note that he uses the notation fx where we use s\. 

For nebular extinction, one often writes the logarithmic extinction as C\ = Cnpil + fx), 
with / H /3 = (e.g., Peimbert & Torres- Peimbert 1977). Then if we adopt s\ — (1 + fx) = 
T \/ r H(3 we have a = Ch/3/0.4343 = r H/ 3. 

For stellar extinction one often adopts the normalization sx = Ax/ Ay = tx/tv (e.g., 
CCM89) in which case a = Ay/1.086 = ry. In order to compare these differently normalized 
nebular and stellar extinctions, one can easily show that 

A J Ay 

Differential extinction between two wavelengths can of course be defined in any system. 
For example, color excesses are differences in magnitudes, like Eb-v — Ab — Ay. These 
can lead to differential extinction curves with new normalizations, like Ex-v/Eb-v — {Ax — 
Ay)/ E B -v — Rv{A-x/A v — 1), accompanied by the introduction of new parameters, like the 
ratio of total to selective extinction, Ry = Ay/Es-y- 



3.1.1. Determining the shape, s x : Stellar extinction 

CCM89 defined the shape of optical stellar extinction analytically using a single param- 
eter, Ry, to describe differences between lines of sight. The shape has been adopted in many 
investigations of the Orion Nebula (Baldwin et al. 1991; Osterbrock et al. 1992, hereafter 
OTV92; Greve et al. 1994), with R v = 5.5. 2 



Optical 

Good optical spectrophotometry of the Orion stars (9 1 and 9 2 ) is presented by Cardelli & 
Clayton (1988) (hereafter CC88) as E X -y/ 'E B _ V . Using measurements from their graphs, 
we show that the optical data (and more scanty infrared data) are fit well by the CCM89 
formula (dash-dotted line in our Fig. 2). The CC88 data (asterisks in Fig. 2) continue to 
about 3 /im" 1 and show that there is a knee in the CCM89 curve at about 2.7 /im" 1 (A3700) 



2 Ry = 3.1 in the diffuse interstellar medium (Savage & Mathis 1979). 
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that is perhaps a bit too high and sharp. Therefore, an analytic modification 3 has been 
made to the CCM89 formula (solid line) to round off the knee smoothly between 2.3 and 
3.3 yum -1 , fitting the CC88 data quite well. 

Ultraviolet 

The ultraviolet data come from IUE measurements by Bohlin & Savage (1981) (hereafter 
BS81) for the same Orion stars (see Fig. 2). The BS81 points (x) do not join on smoothly to 
the CC88 optical data. CC88 reanalyzed the ultraviolet data, showing how it could reason- 
ably be joined smoothly to their new accurate optical spectrophotometry (asterisks). The 
main effect on the BS81 data (x) is to move it up slightly (~ 0.25 mag shift), concomitantly 
increasing the expected ultraviolet extinction (see Fig. 2). 

Note that in the ultraviolet, the analytic CCM89 formula does not fit either the CC88 
or BS81 observations: the predicted bump is too strong. However, it is possible to adjust a 
few constants 4 in the CCM89 formula to obtain a modified CCM89 curve (solid line) that 
runs smoothly through the CC88 data. 

To gauge the importance of these modifications for the determination of extinction- 
corrected fluxes, consider the situation for the N n] A A2 140-43 pair near the peak of the 
extinction curve. The unmodified CCM89 curve exceeds ours by SE\_ v /Eb-v = 0.85. For 
a typical value Eb-v — 0.3 the difference would be St = 0.23; thus the N il] line would be 
26% (0.1 dex) stronger if corrected according to the unmodified CCM89 curve. 

Near infrared 

Near 1 /xm the CCM89 A\/Ay curve becomes a power law in A -1 of index 1.6. Most stellar 
curves exhibit this "common power law" region of the extinction curve, independent of Ry 
(Martin & Whittet 1990). The latter analysis suggests an index of 1.85, which we adopt in 
our extinction curve. This modification will not have an impact here, as the observations 
analyzed do not extend much beyond 1 /im. 

We have included a table of fx values for a number of observed emission line wavelengths 
(see Table 2), calculated directly from our modified set of CCM89 equations (with R v = 5.5). 



3 CCM89's equations 3a and 3b are modified to a(x) = 1.192528 + 0.27592(a; - 2.3) - 0.15733(a; - 2.3) 2 
and b(x) = 0, but only for 2.3 [ivcT 1 < x < 3.3 ixmT 1 . 

4 CCM89's equation 4b becomes: b(x) = -2.9 + 1.825a; + 0.93/[(a;-4.65) 2 + 0.263] + F b (x) for 3.3 ^m" 1 < 
x < 8 ^m -1 . 
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For other wavelengths, refer to CCM89 and our analytic modifications discussed above. 

3.1.2. Setting the amplitude, a: Optimal fitting of emission lines 

Measurements of emission lines in nebulae are often used to gauge the extinction (e.g., 
CP70; OTV92; Greve et al. 1994; Esteban et al. 1998, hereafter EPTE98; Esteban et al. 2004, 
hereafter EPG04). If the shape s x of the extinction curve is known, then the amplitude a 
in equation 2 can be determined with even a single line pair from the same ion, assuming 
the relevant emissivities j\ are known. The Balmer and Paschen lines are popular as their 
emissivities are well described by H I recombination theory and can be computed readily 
for various electron temperatures (T e ) and densities (N e ). See Storey & Hummer (1995) for 
detailed access to the results of their computations using an interactive data server. For the 
illustrations below we adopt case B emissivities with T e = 8000 K, N e = 2500 cm -3 for FOS- 
1SW and STIS-SLITlc (as determined from nebular models presented in § 4.4), T e = 9000 K, 
N e = 5000 cm- 3 for the remainder of the STIS slits, and T e = 8300 K, iV e = 8900 cm" 3 
for the EPG04 observations (as determined from their temperature and density diagnostic 
ratios). Lines from other ions can be useful too (e.g., He I, [O n], [S il]), as further illustrated 
in § 4 (He i) and § 3.3 ([O n], [S n], He i). 

Method 

Sometimes the amplitude is taken from only one line ratio, like Ha/H/3, but often other lines 
are available too. Past methodology has been to form various line ratios, compute the ampli- 
tude (e.g., the equivalent A v ) from each independently (sometimes with discrepant results), 
and then compute some average. A new approach adopted here is to fit the observed line 
fluxes directly to equation 2, not forming line ratios at all. The fit is carried out straightfor- 
wardly by non-linear least squares, with two unknowns: a and a constant multiplier (which 
is proportional to the ionic abundance through the emission measure.) In addition, from 
the goodness of fit of the data to the model, we obtain the la confidence intervals on the 
parameters. This approach has some advantages: avoiding the inevitable bias from errors 
in the chosen normalizing line(s), using a standard treatment of the measurement errors 
to weight the fit and hence avoiding the ambiguity associated with deciding what average 
amplitude to compute. 

If later one wants to tabulate extinction-corrected line ratios relative to, say, H/5 (or 
He I 4471 if He I lines are used in the extinction fitting), then the extinction model gives the 
best estimate of the corrected flux for this normalization, optimally consistent with all other 
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lines in the series. When the He I recombination lines are fit simultaneously, then the optimal 
He + /H + ratio is obtained, along with its formal error, without reference to particular He I 
to H I dereddened line ratios. 



3.2. Comparison between stellar and nebular extinction curves 

The CP70 nebular extinction curve is an alternative commonly used in the calculation 
of the extinction correction (EPTE98; EPG04). This curve was determined from optical and 
radio observations of four lines-of-sight in the Orion Nebula. The radio (1.95 cm) observations 
are used to determine the absolute extinction, A Hl g, and the piecewise normalized extinction 
curve, 1 + fx, is obtained for various lines: 

h = (4) 

Calculations of C^p offer a simple means by which to compare our modified CCM89 
curve with that of CP70. Performing a least-squares analysis on EPG04 data using the 
CP70 curve (and interpolations thereof), we get Ch/? = 0.65 ± 0.03, as compared to Ch/3 = 
0.85 ±0.04 when we use our extinction curve. EPG04 had C H/ 3 = 0.76 ±0.08 when using this 
same series (up to H15, A3712) of H I lines. However, this comparison of C H/ 3 is not enough. 
The shape of the CP70 curve is investigated more thoroughly in the following sections. 

3.2.1. Infrared discrepancy 

The zero-point is calculated by CP70 using radio data at (more or less) the same spatial 
position as the optical data. The rest of CP70's fx is determined directly from specific 
Balmer or Paschen lines (as tabulated in Peimbert & Torres- Peimbert 1977), so that most of 
the curve plotted in Figure 3 is interpolated (by CP70; Peimbert & Torres-Peimbert 1977; 
Walter 1993). 

Near 1 /xm -1 , the CP70 curve is dramatically below our extinction curve (see Fig. 3). 
This is getting into the "common power law" region of the extinction curve mentioned in 
§ 3.1.1. Nowhere is there evidence for the strong curvature implied in the CP70 extinction 
curve in the region from 1 /im" 1 to zero wavenumber. The IR discrepancy might simply 
be the result of an inappropriate interpolation of the CP70 extinction curve. Looking at 
Figure 1 of CP70, it appears as if other interpolations to the V band are possible. If the 
interpolated fy were a bit more negative (—0.12 or —0.13 instead of —0.094 would seem 
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reasonable) and the CP70 curve (shown in our Fig. 3) was refit (and rescaled, see below) , it 
would be noticeably closer to the stellar one from the optical through the near infrared. 

The change in the interpolated f v would affect R v as derived from 

Rv = 4±^, (5) 

JB - JV 

resulting in R v = 4.5 (f v = -0.12) or 4.2 (f v = -0.13) compared to the CCM89 and 
CP70 value of Ry = 5.5. However, the fx depend on comparing radio and optical data (see 
equation 4). This radio/optical comparison is almost certainly affected by a 10-20% reflection 
effect in the optical (e.g., Henney 1998). Because of reflection, the lines like H/3 are measured 
to be too bright relative to the radio, resulting in a low estimate of the absolute extinction, 
Ah/3- In the optical where the differential extinction Ex-up is not affected much by reflection, 
fx = Ex-ap/ Aup will be too high. A typical logarithmic extinction Cup (Ch/3 = 0.4A H /?) is 
0.5, so that with f v = —0.12, the logarithmic extinction at V is 0.5 x (1 — 0.12) = 0.44 or 
Ay ~ 1. A 20% reflection effect at V would cancel out real extinction at the level A Ay = 0.2, 
lowering what should be 1 to 0.8 (apparent). Then 0.8 x 5.5 predicts Ry (apparent) = 4.4. 
Both the sign and order of magnitude of the reflection correction appear to be right. Note 
that there are other possibilities of systematic error as the radio/optical comparison also 
depends on the adopted temperature, matching the spots observed, etc. 

To remedy this artificial alteration of Ry, we introduce a correction to A-^p - effectively 
a rescaling of fx - to maintain Ry = 5.5. With fy = —0.12 in our new interpolation of the 
CP data, one needs to divide all fx by a factor 1.198 to ensure that Ry = 5.5. This rescaled 
CP70 absolute extinction curve is presented in Figure 3. 

The shape of our rescaled CP70 curve in the visible is similar to that of our modified 
CCM89 curve and so regardless of which of these is used, one should be able to get good 
differential fits to a series of optical lines. The differential extinction between Ha (1.52 jum" 1 ) 
and H<5 (2.44 /im -1 ) is the same in both curves, and one obtains more or less the same Ch/3 
when fit to the Balmer lines, independent of the choice of curve (see § 3.3.1). 



3.2.2. Ultraviolet 

Walter (1993) extended the CP70 curve to the ultraviolet by attaching on the BS81 
data assuming Ry = 5.5. This extended extinction curve was also described by Rubin et al. 
(1993). We have derived fx values from the original BS81 E X -v I 'Eb-v form (as in Rubin 
et al. 1993) for both the CP70 and rescaled CP70 curves, generating a shape of the ultraviolet 
extinction for each, as seen in Figure 3. This extension lies below the CC88 revision of the 
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BS81 data in the ultraviolet and does not join smoothly to the optical portion of the CP70 
curve. One of the reasons is that the CP70 curve has been extrapolated to A3500 (see their 
Table 2), incorrectly high, beyond the highest Balmer transitions that were observed (H9, 
A3835). When combined with the ultraviolet data (BS81 version) this really sticks out as 
a kink in the curve at A3500; a smooth join, like all (stellar) extinction curves observed, is 
not possible. The reddening correction of a line like [O n] A3727 would be in doubt. Our 
modified CCM89 curve seems a more reliable alternative in this "extrapolated" region. 

3.3. The stellar extinction curve: validation 

It is important to address whether the shape of extinction derived from stars is ap- 
propriate for diffuse nebular emission, where the presence of scattered light might be an 
issue. We have therefore carried out extensive tests of this technique using consistent sets of 
nebular measurements of hydrogen lines in Orion, both from the literature (CP70; Peimbert 
k Torres-Peimbert 1977; OTV92; EPTE98; EPG04) and from our HST and ground-based 
(Baldwin et al. 2000, hereafter BVV00; Blagrave et al. 2006a,b) observations. 

Here we exploit our HST FOS and STIS (SLITlc) observations for the line of sight 1SW 
(see Table 1) as, with the addition of several UV lines, they offer an extensive spectral range 
for the validation. As will be discussed here, and again in § 4.6 (Figures 9 and 10), the data 
are well fit using our extinction curve. 

3.3.1. Optical 

In addition to our own STIS and FOS observations we have used others' observations 
to examine the fit of the stellar-derived extinction curve in the optical. Using the extensive 
observations of EPG04, we have determined the best-fit extinction curve (see § 3.1.1) using 
the unblended H I Balmer lines (up to H17, A3697) and Paschen lines (up to P18, A8438), 
resulting in Cu/3 = 0.82 ± 0.04. The resulting curve and residuals from this fit are shown 
in Figure 4. Over the spectral range that can be accessed by these lines, the fit is good, 
demonstrating that our extinction curve provides a good empirically-derived interpolation 
formula for differential extinction. 

The EPG04 data were also fit using the scaled form of CP70 developed in § 3.2 and 
the same H I lines as above. The resulting curve and residuals are included in the bottom 
panel of Figure 4. We find Ch/3 = 0.76 ± 0.04, in agreement with Ch/3 = 0.76 ± 0.08 found 
by EPG04. As shown in the top and bottom panels of Figure 4 and discussed in § 3.2, 
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the differential extinction is well represented by both our scaled CP70 and modified CCM89 
curves. 

We have also carried out fits of the various data sets using our curve and R v = 3.1, 
appropriate to the diffuse interstellar medium (Savage & Mathis 1979). By comparison to 
the Orion curve (with Ry = 5.5) this curve is steeper throughout the optical, producing 
relatively less (more) extinction compared to H/3 in the near- infrared (blue). This produces 
a markedly inferior fit. 

3.3.2. Ultraviolet 

Validation using nebular observations is especially important for the ultraviolet portion 
of the extinction model. The STIS and FOS spectra include the common upper level pair 
[O n] AA2471, 7325 which can be used in this validation. 

The [O n] A2471 line is actually a blend of transitions from 2 P\/2 and 2 Pz/2 to 4 S , 3 / 2 
(Osterbrock 1989). Transitions from these upper levels also lead to the formation of a pair 
of blended near- infrared lines. Each upper level contributes a line to a blend at A7321 (tran- 
sitions to 2 D 5 / 2 ) and another line to a blend at A7332 (transitions to 2 D 3 / 2 )] the combined 
near-infrared line will be denoted A7325. Thus the 12471/^7325 line ratio is not quite the ideal 
case with a single common upper level in which the emissivities of the line pair are simply 
proportional to the radiative transition probabilities divided by the respective wavelengths. 
Instead, this line ratio depends on the collisional excitation of 2 Pi/2 relative to 2 Ps/2- How- 
ever, calculations with a five-level atomic model show that the line ratio is not very sensitive 
to T e and N e over the range expected in the Orion Nebula. We adopt 0.75 using modern 
collision strengths (Pradhan et al. 2006) and older transition probabilities (Zeippen 1982). 
We compare these results with a second ratio, 0.81 derived using the modern set of transition 
probabilities recommended by Wiese et al. (1996). These are both calculated for T e = 8000 K 
and N e = 2500 cm" 3 (see § 4.4). 

In Table 3, the predicted /2471/J7325 ratios are compared to our observed ratio corrected 
using our extinction curve and H I lines. Use of the modern transition probabilities results 
in a 15-20% over-prediction, whereas the older transition probabilities result in a slightly 
lower 10-15% over-prediction. These modern transition probabilities have been found to 
yield anomalous density and temperatures from [O 11] lines (see Wang et al. 2004; Blagrave 
et al. 2006a, EPG04) and thus it is not unexpected that the older transition probabilities 
yield a slightly more consistent result. 

As a check on the theoretical predictions, we have also examined the [O 11] /7321/ ^7332 



- 12 - 



line ratio which is slightly more sensitive to the relative upper level populations, but not 
sensitive to the applied reddening correction. The prediction is 1.25 using Zeippen (1982) 
transition probabilities and 1.23 using Wiese et al. (1996) transition probabilities. Our ratios 
are 1.21-1.24 for the STIS observations. Other observations in Orion by OTV92 give 1.134, 
while BVVOO find 1.25 and EPG04 find 1.30. Our ability to accurately predict (within 15%) 
both the /7321//7332 line ratio and the common- upper-level line ratio /2471 / '^7325 supports the 
validity of our analytic extinction curve in the UV. However, it would be beneficial to have 
further constraints in the UV. This is addressed in § 4. 

3.3.3. Near infrared 

Similarly, the [S 11] common upper-level line pair Jio3oo/-^4072 can be used to confirm the 
validity of the extinction curve in the infrared. Our data do not extend into the infrared, 
but those of EPG04 do. The three observed lines at A 10300 and two at A4072 have been 
corrected for extinction using our best-fit extinction curve (see § 3.3.1). Since this is a ratio 
being used for validation and is not in the fit, the visible line (which tends to have the lower 
error) has been set to I corrected/ 1 predicted = 1 in Figure 4. The agreement is good. 

As pointed out in Porter et al. (2006), EPG04 have also observed a series of He I common 
upper-level line pairs. These line pairs have been overplotted on Figure 4 too. The IR pair 
members have I corrected/ {predicted = 1-2 ± 0.5. While there are large differences between the 
observed and predicted line strengths, there is no significant systematic bias, indicating that 
our extinction curve is also valid between the near-IR and near-UV. 



4. Anomalous He I decrements 

Based on predictions from recombination theory we have used the H I Balmer and 
Paschen series to calibrate the reddening curve. The same should be possible using recom- 
bination theory (Smits 1996; Benjamin et al. 1999; Porter et al. 2005) and observations of 
the He I lines (see energy- level diagrams in Fig. 5). 

4.1. Extinction-corrected data 

Table 4 presents our He I data from FOS-1SW and (the section of our STIS spectrum 
that covers 1SW) STIS-SLITlc observations as well as complementary data from compre- 
hensive ground-based studies by OTV92, EPTE98, BVVOO, and EPG04. 
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Each of the He I lines has been corrected for reddening using our best-fit curve as 
normalized using the observed H I Balmer (and Paschen, where available) series lines, and 
each line's case B prediction. Case B predictions depend on T e and N e (Storey & Hummer 
1995) and so are unique to each observation set: T e = 9000 K, N e = 4000 cm -3 (OTV92); 
T e = 9500 K, iV e = 5700 cm" 3 (EPTE98, position 2); T e = 9000 K, N e = 5000 cm" 3 
(BVV00); T e = 8300 K, N e = 8900 cm" 3 (EPG04); T e = 8000 K, N e = 2500 cm" 3 (FOS- 
1SW, STIS-SLITlc). 

As is common, the line strengths are given relative to the singlet line A4471. 

The quoted uncertainties for OTV92 and EPTE98 are estimated from their description 
of the quality of the data. Those of EPG04 are directly from the per cent errors in their 
Table 2. BVV00 and FOS-1SW uncertainties were determined using the signal to noise ratio, 
S/N, as found from fitting each emission line to a five-parameter (area, central wavelength, 
FWHM, continuum baseline and continuum slope) single Gaussian model. This S/N is 
combined with the integrated line flux to produce the uncertainty. The STIS data is similarly 
fit with a model which takes into account the 0'/5 STIS slit width (see Appendix A in Rubin 
et al. 2003). Again, the resultant S/N from this fit is combined with the integrated flux to 
determine the uncertainty. 

4.2. Case B 

He I case B predictions for FOS-1SW and STIS-SLITlc T e and N e (n < 5, Porter et al. 
2005; n > 5, Smits 1996) are shown in column (3) of Table 4 and are quite representative 
for the range of nebular densities and temperatures we consider. 

Most of the reddening-corrected singlet and triplet lines (cols. (6)-(ll)) are well-approximated 
by case B recombination alone (col. (3)). 

In comparing their case B results with the OTV92 observations, Benjamin et al. (1999) 
noted (their Figs. 2a and 2b) that the lines (singlets and triplets) shortward of A3820 were 
consistently brighter than predicted, increasingly with decreasing wavelength, except A3 188 
which appeared to be accurately predicted. They suggested that this upward trend shortward 
of A3820 could be due to unaccounted for radiative transfer effects, though this seems very 
unlikely given the lines involved. Inclusion of radiative transfer (see discussion below) would 
cause the plotted A3188 ratios to increase - possibly enough to agree with the apparent 
upward trend shortward of A3820 - while not affecting the other lines shortward of A3820. 
If alternatively this upward trend arose from inadequate/unmodeled extinction corrections, 
there would have to be a very sudden (relative) decrease in the extinction in order for 
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the extinction-corrected lines not to appear so bright. It is informative to note that the 
discrepancy in this instance cannot be explained by the inaccurate extrapolation of CP70 
as discussed in § 3.2.2 (since CP70 is not used in the extinction correction). It seems 
most probable that there is a systematic problem with the absolute calibration used and/or 
significant measurement errors due to strong atmospheric extinction varying rapidly with 
wavelength for A < 3500 A (OTV92). See also § 4.5, Figure 8. 

In their discussion of nebular model predictions of the He I lines compared to the EPG04 
observations, Porter et al. (2006) have noted that He I lines of the 2 1 S — n l P series in 
the same wavelength range (AA3448, 3355, 3297) are also brighter than the expected case B 
predictions. Analyzing the line ratios between successive members of this series (their Figs. 1 
and 5) they conclude that there are problems with these observational results, possibly with 
the extinction (over-)correction. The same observations, corrected according to our best-fit 
extinction curve, are presented in Table 4, where it can be seen that the observations and 
case B (or our full models) are in agreement within the errors. 

As seen in Table 4, use of our extinction curve brings A7281 into agreement with the 
case B predictions (as suggested by Porter et al. 2006). Many of the other near-IR lines ob- 
served by EPG04 are still scattered about the best-fit expectation (Table 4 and graphically in 
§ 3.3.3, Figure 4). These infrared lines are difficult to observe because of atmospheric effects 
(OTV92) and the errors might be underestimated (Porter et al. 2006), but the important 
thing to note is that there does not appear to be any systematic discrepancy - except as 
noted below for the infrared lines of the triplet series 3 3 S — n 3 P. 

4.3. Anomalous decrements 

Despite the validity of case B for the prediction of most lines, there remains a systematic 
mismatch between observed and predicted ratios for a number of the He I lines from the series 
2 3 S - n 3 P, 2 3 P - n 3 S, 3 3 S - n 3 P, and 2 1 S - n 1 P. 

For the triplet lines, the underlying reason is the metastability of the 2 3 S term which 
can only be depopulated via photoionization, through collisional transitions to 2 1 S' and 2 : P, 
or through the (strongly forbidden) radiative transition 1 1 S — 2 3 S (Osterbrock 1989). The 
visible/near-UV lines of the 2 3 S - n 3 P series (including AA3889, 3188, 2945, 2829, etc.) 
are affected by self-absorption from the metastable 2 3 S term to the corresponding n 3 P 
term. This results in case B recombination theory over-predicting these lines (see Table 4). 
Note that although A3889 is blended with a Balmer line, it is still possible to predict the 
observed intensity of the He I component and use it as a valuable diagnostic. For the entries 
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in Table 4 we used case B H I predictions and the He + /H + derived from the other lines to 
deblend A3889; this also enables a "case B" prediction for the blend (see col. (3) in Table 4). 
The space-based observations that we report here are particularly valuable for verifying the 
expected diminution of this self-absorption effect in higher members of the series. This 
appears to be the case as discussed in the modeling below. 

For self-consistency, there must also be a detectable resonance fluorescence enhancement 
of the 2 3 P — n 3 S series as some of the electrons promoted by self-absorption to n 3 P 
radiatively cascade back down to 2 3 S' via alternate routes. Indeed, this is quite evident 
in the 2 3 P — n 3 S series, primarily in A7065 (n = 3), but detectable in A4713 (n = 4) 
and A4121 (n = 5) too (see Table 4). Similarly, we expect to see an enhancement of other 
triplet series including the near-IR series, 3 3 S — n 3 P. The data in Table 4 show this 
enhancement. We will return to analysis of these series in § 4.5. 

The singlet lines of the series 2 1 S — n 1 P, A5016 (n = 3) (and to a much lesser 
extent A3965, n — 4) are also observed to be consistently less than the case B prediction. As 
discussed in Porter et al. (2006), this is most probably a result of UV lines (A537, 1 1 S — 3 1 P 
and A522, 1 1 S' — 4 1 P) escaping the nebula, resulting in case B being shifted slightly toward 
case A (as defined in Baker & Menzel 1938). 

4.4. Photoionization models 

The CLOUDY photoionization code (Ferland et al. 1998) accounts for many radiative 
transfer effects, including those relating to the metastable 2 3 S term. It is also possible to 
run CLOUDY models that exclude all line transfer in order to investigate its importance. 
With this in mind, a series of CLOUDY photoionization models was developed. 

Self-absorption effects depend on the line width because this affects the opacity. The 
He I lines (among others) observed in the ground-based echelle spectra (BVV00; Blagrave 
et al. 2006a,b) have line widths (FWHM) in excess of the Doppler/thermal and instrumental 
widths. When the Doppler/thermal and instrumental widths are subtracted (in quadrature) 
from the measured FWHM there remains an unaccounted for contribution to the broadening 
which has been labeled as "turbulence". In the case of the He I singlet lines (which do not 
have any fine-structure levels and therefore will not have any further FWHM contribution 
from a line blend) in the 1SW ground-based echelle spectra (Blagrave et al. 2006b), we find 
FWBM turb = 15.5 ±2.4 km s -1 . This is quite consistent with the "turbulence" contribution 
quoted in O'Dell et al. (2003) for the He I A5876 triplet line: 18.4 ± 2.9 km s _1 . This latter 
line and other triplet lines have another broadening contribution brought on by transitions 
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in the fine structure (J levels) and therefore are slightly less reliable in the determination of 
the additional "turbulence" contribution. As was mentioned in O'Dell et al. (2003) and is 
also found here, the magnitude of the "turbulence" contribution to the line width appears 
to decrease as one probes along the line of sight from cooler (associated with H I and He i) 
to hotter (associated with [O ill], [O n], [S n], etc.) regions of the nebula. As the He I 
recombination lines preferentially probe the cooler region of the nebula, the "turbulence" 
that these lines reveal is an upper limit for the nebular model. 

This turbulent contribution to the line width will affect the magnitude of the radiative 
transfer effect on the 2 3 S — n 3 P series of lines (and the related lines, like A7065), as it 
broadens the lines and lowers the effective optical depth and line trapping. The effects 
from radiative transfer should diminish as the turbulence is increased. To investigate the 
effects of this change in turbulence on line predictions, we developed a series of constant 
density (10 4 cm" 3 ) and constant temperature (10 4 K) CLOUDY models each with a different 
turbulence parameter, v tur b = FWHM teri) / V4 In 2 (as defined in CLOUDY), and excluding 
"induced processes" (as defined in CLOUDY, these include continuum fluorescent excitation 
and induced recombination). Then each of these constant-turbulence models was varied in 
thickness, in order to develop plots of line ratios as a function of To(A3889) - the optical depth 
at the center of He I A3889 for zero turbulence velocity (see Fig. 6). For the normalizing line 
we use A4471 which is not significantly affected by radiative transfer. Our Figure 6 can be 
compared with Figure 4.5 of Osterbrock (1989), which shows similar variation, but in the 
latter case as a result of changes in the nebular expansion velocity. At t = all models 
reduce, as is expected, to their case B values. 

Figure 7 plots the line strengths, h/hun, as a function of each model's "turbulence" pa- 
rameter for fixed r (A3889)~ 27 which is appropriate to the column density of our full nebular 
model below. In this case to enable the best comparison between the CLOUDY model predic- 
tions and our observations, the "induced processes" are included (regardless, in our full model 
the relative importance of the induced processes is much diminished because the nebula be- 
comes optically thick to the continuum radiation responsible for these excitations). From 
the overplotted observed /a7065 / Ixaah in Figure 7, we would predict v tur b = 20 — 30 km s _1 , 
whereas from A3889 (deblended), A3188 and A2945, v turb = 10 — 12 km s _1 . This latter 
prediction is consistent with the "turbulence" contribution (FWHM tur6 ~ 16 km s _1 ; i.e., 
Vturb ~ 10 km s^ 1 ) which was determined from the observed He I lines' FWHM above. The 
value vturb = 12 km s _1 will be adopted in the CLOUDY models discussed hereafter. Models 
with this value of the turbulence should be able to roughly explain the observations of the 
hitherto anomalous He I lines. As expected, the results of the full nebular models to be 
discussed (Table 4) are slightly different than the results for constant N e and T e in Figure 7, 
in particular giving even better agreement for A7065. 
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Our full nebular models are similar to the Orion Nebula model discussed in Baldwin 
et al. (1991) - i.e., a closed geometry and constant pressure. Two such models (see Tables 5 
and 6) were created, as was done in Blagrave et al. (2006a): one with a Mihalas stellar atmo- 
sphere model (Mihalas 1972) and H II region abundances as defined in CLOUDY (Baldwin 
et al. 1991; Rubin et al. 1991, OTV92) (model M); and one with a Kurucz stellar atmo- 
sphere model (Kurucz 1979) and EPG04 Orion Nebula abundances (model K). The Mihalas 
(non-LTE) stellar atmosphere has been shown to best represent the incident continuum 
radiation (BVV00), but to test the robustness of our results we also include the Kurucz 
(LTE, line-blanketed) stellar atmosphere in our second model. Each of these models has 
Vturb = 12 km s -1 . The other parameters (density, radius, etc.) were varied so as to best re- 
produce the SLITlc H/9 surface brightness and the temperature- and density-sensitive lines, 
as summarized in Table 7 for both models. 

The models' predictions of the He I lines given in columns (4) and (5) in Table 4 include 
the deviation from case B introduced by the metastable 2 3 S term. The effect of this 
metastability is reflected in the models' (non-zero) line optical depth, r(A3889) = 13. Note 
that He I A3889 (n = 3) appears as a blend with H I A3889 and He I A2829 (n = 6) appears as 
a blend with [Fe iv] A2829. Neither of these blends is included in our quantitative analysis, 
but the predicted value for each is shown in a separate row in Table 4. 

Models M and K were modified to exclude radiative transfer in order to further inves- 
tigate the validity of the CLOUDY models. This exclusion has the expected effect on the 
lines: returning the triplet line predictions to roughly their case B values and those singlets 
connected radiatively to the ground state are reduced to their case A values. 

4.5. Comparison with observations 

We have modeled the extinction correction for OTV92, EPTE98, BVV00, EPG04 and 
our FOS-1SW and STIS-SLITlc observations using our extinction curve and the H I Balmer 
(and Paschen, where available) lines and a subset of He I lines, excluding lines from the 
triplet series 2 3 S - n 3 P (AA3889, 3188, 2945, 2829), 2 3 P - n 3 S (AA7065, 4713, 4121) 
and 3 3 S — n 3 P (all lines) and the singlet line, A5016. Our Cup from the fit of H I 
and the subset of He I lines is consistent with that found using H I lines alone (see end of 
Table 4), suggesting that both H I and this subset of He I lines are reliable determinants of 
the extinction correction. 

For comparison with these extinction-corrected observations, I C orrected, we adopt sepa- 
rately predictions, I pred icted, from case B and CLOUDY model M. Values of \og(I corrected / I predicted ) 
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are plotted in Figure 8 for all observed members of four series for both case B and model M 
predictions. Many data sets are shown, but EPTE98 is excluded simply to avoid visual 
clutter. For the 2 3 P — n 3 D series, which includes A4471 (n = 4), case B and model M 
are in close agreement. This triplet series is not significantly affected by radiative transfer 
effects and the theory and observations are in close agreement. The other triplet series are 
clearly affected, by self- absorption (2 3 S' — n 3 P) or by resonance fluorescence enhancement 
(2 3 P — n 3 S, 3 3 S — n 3 P). Case B is a poor approximation, but the agreement between 
theory and observations improves dramatically as we switch from case B to model M. (x 2 is 
reduced by a factor of roughly 10.) 

Two features deserve further discussion. The OTV92 data shortward of 3700 A (includ- 
ing A3188) most likely suffer from systematic errors (§ 4.2). This can be seen graphically in 
the 2 3 P — n 3 D series of Figure 8 when comparing OTV92 (circles, no uncertainty bars) 
and others' datasets; the OTV92 data points (with one exception, A3479) all lie significantly 
apart from the other data. Second, A9464 measured by EPG04 appears to be too faint. This 
line has a number of neighboring near-IR H I Paschen lines and He I lines which appear to 
be well-fit by our extinction correction. There would have to be a large sudden increase in 
the amount of extinction for the extinction correction to be able to account for this anomaly. 
This suggests that A9464 may have a larger than reported measurement uncertainty (see 
also § 4.2 and Porter et al. 2006). 



4.6. Extinction corrections revisited 

The STIS and FOS data are corrected for extinction using our extinction curve, the 
unblended H I lines, the subset of He I lines and the common-upper-level pair of [O n] lines 
(refer to § 4.5 and § 3.3.2). These are presented in Table 8 relative to the Hf3 predicted from 
the fit, along with the derived Ch/3- 

In Figure 9, we show the extinction-corrected FOS-1SW line intensities compared to 
the case B theory (the residuals of the fit) as well as the extinction correction applied 
(here expressed differentially with respect to Cup)- The He I lines AA7065, 3188, and 2945 
are excluded from the fit, but their values relative to case B and model M predictions are 
overplotted. Model M provides a better prediction for the lines affected by radiative transfer. 

The analysis for the SLITlc data is shown in a similar way in Figure 10. The He I lines 
AA7065, 5016 and 2945 (and H I line A3704, because of large uncertainty) are excluded from 
the fit, but their values relative to case B and model M are overplotted. Again, model M 
provides a better prediction for the lines affected by radiative transfer. 
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In fitting the UV extinction curve, we conclude that A3889 (a blend in any case), A3188, 
A2945, and A2829 (also a blend) should be excluded - pending correction for radiative transfer 
effects - but that lines associated with higher n (i.e., n > 6) can be included to ensure a 
more robust determination. The visible/IR extinction correction can be strengthened by 
including case B predictions for He I triplet and singlet lines, but as discussed above, we 
conclude that lines from the series 3 3 S — n 3 P (all lines), 2 3 P — n 3 S (AA7065, 4713, 
4121) and 2 1 P — n 1 S (A5016) should be excluded - again pending correction for radiative 
transfer effects. 

4.7. He+/H+ ratio 

In fitting the H I and He I lines simultaneously, we are also able to calculate the He + /H + 
abundance directly from the parameterization used in the fitting (see § 3.1). These abun- 
dances are included at the bottom of Tables 4 and 8. Table 4 also shows the abundances 
from the original reference (when published). Our new analysis agrees with these within the 
errors. Our confidence intervals reflect implicitly the consistency of the extinction correction 
as well as the match to theory and observational errors and so are somewhat larger. 

The He + /H + ratio could in principle vary with position in the nebula if the ionization 
structure (hence the ionization correction factor, ICF, used to convert He + /H + to He/H) 
differs on account of nebular structure. However, this is likely to be a small effect. It is the 
case that the He + /H + values for the different positions agree with one another within their 
errors. The unweighted mean and its standard deviation is 0.0883 ± 0.0019. If weighted, the 
ratio is 0.0882 ± 0.0004. Neither accounts for any systematic error. 

Extensive discussion of the ICF required to determine the atomic He/H ratio would take 
us beyond the scope of this paper. However, our CLOUDY model M gives 1.14 and model K 
gives 1.00, which are within the range given by other authors; 1.1 is a typical value with an 
uncertainty of roughly ±0.05. Adopting this we obtain He/H = 0.097 ± 0.002 ± 0.005. The 
latter uncertainty from the ICF is probably the major source of uncertainty, greater than 
that in measuring the ionic abundance. 

5. Summary 

Using slight modifications to the valuable stellar extinction curve developed by CCM89, 
we have been able to develop an accurate new analytic method of determining the nebular 
extinction curve over an extensive wavelength range in the Orion Nebula. This curve has 
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been rigorously tested with currently available near-IR, optical, and ultraviolet ground- 
based and space-based data, standing up as a robust measure of the extinction. We have 
also compared this new curve with the CP70 nebular extinction curve. We have confirmed 
that the discrepancy, with respect to theoretical expectations, in some near-IR and near-UV 
He I lines measured by EPG04 is a result of inaccurate extinction correction, but that the 
UV discrepancy in observations by OTV92 is a result of calibration and/or measurement 
uncertainty and not an extinction (or radiative transfer) effect. 

On the foundation of this new extinction analysis, we have measured systematic anoma- 
lous He I decrements, compared to case B, associated with the 2 3 S — n 3 P, 2 3 P — n 3 S and 
3 3 S — n 3 P series, and to a lesser extent with the 2 1 S — n 1 P series, none of which can be 
explained by any adjustment to the extinction curve. Qualitatively, these anomalies are as 
expected from radiative transfer effects, for the triplets arising from the metastability of 2 3 S' 
(Osterbrock 1989). Furthermore, modeling of the radiative transfer effects using CLOUDY 
produces a remarkable quantitative agreement between theory and observation. 

Because He I case B recombination theory is not reliable for a subset of He I permitted 
lines associated with these aforementioned series, those lines most affected (A3889, A3188, 
A2945, A2829, A7065, A4713, A4121, A5016 and those of the triplet series 3 3 S - n 3 P) 
must be excluded in the determination of the amplitude of the extinction curve, or in the 
calculation of He + /H + abundance. Alternatively, they could be included after adjustment 
from their case B predictions by modeling radiative transfer effects. 
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Table 1. HST STIS V2 and FOS observations prior to reddening correction I^p) 
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o 
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o 
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II I 
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0011 
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0005 
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0276 
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H I 
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o 


0361 


± 


o 


0010 


o 
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3889 
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He I+H I+[Ne III] 
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4026 
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4170 
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3777 
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0091 
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0002 





0074 


± 





0003 
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He I 


4388 





0048 
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0002 





0045 


± 





0001 





0039 
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0048 


± 
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0042 ± 0.0006 


He I 


4471 
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0005 
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0005 
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0416 zb 0.0008 


H I 


4861 


1 


0000 


± 
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1 


0000 


± 





0110 


1 


0000 


± 





0073 


1 


0000 


± 





0131 


1 


0000 ± 0.0081 


He I 


4922 





0093 


± 





0010 
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± 





0002 





0142 


± 





0010 





0106 


± 





0007 





0121 ± 0.0013 


[O HI] 
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1 


2387 


± 
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1 
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± 
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1 
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± 





0096 


1 
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± 
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[O m] 
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5518 





0049 


± 





0003 





0039 


± 





0001 





0028 


± 





0002 





0036 
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± 
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He I 
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± 





0020 
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± 





0018 





1643 


± 





0010 
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± 
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1449 zt 0.0013 


[N II] 


6548 





2341 


± 





0031 





2667 


± 





0031 





2101 


± 





0021 





2461 


± 





0035 






H I 


6563 


4 


4094 


± 





0510 


3 


9699 


± 





0498 


3 


9989 


± 





0294 


4 


4652 


± 





0665 


4 


2067 ± 0.0336 a 


H I 


6563 










































3 


7418 ± 0.0265 b 


[N II] 


6583 





7209 


± 





0102 





8343 


± 





0098 





6460 


± 





0071 





7654 


± 





0106 






He I 


6678 





0549 


± 





0008 





0504 


± 





0007 





0476 


± 
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0539 


± 





0009 





0439 ± 0.0031 a 


He I 


6678 













































0433 zt 0.0016 b 


[S ii] 
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± 





0007 
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± 
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± 
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± 
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[S ii] 
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± 





0012 





0636 


± 





0008 
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± 
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± 
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He I 


7065 
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± 
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± 
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± 
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± 
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2457 


± 





0029 





2298 


± 





0020 





2742 


± 
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C ii 


7236 





0048 


± 





0004 





0041 


± 





0002 





0038 


± 





0002 





0041 


± 





0004 






O I 


7254 





0022 


± 





0002 





0019 


± 





0001 





0014 


± 





0001 





0022 


± 





0001 






He I 


7281 





0119 


± 





0003 





0112 


± 





0001 





0095 


± 





0002 





0113 


± 





0002 






[O II] 


7319 





1144 


± 





0014 





1547 


± 





0019 





0896 


± 





0010 





1215 


± 





0020 





3344 ± 0.0056 


[O II] 


7330 





0936 


± 





0012 





1276 


± 





0015 





0723 


± 





0008 





1006 


± 





0016 




Blend 



a FOS/G780 grating 
b FOS/G570 grating 
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Table 2. Orion extinction curve presented as f\ for select emission line wavelengths 



A (A) x (urn" 1 ) ID f\ 



1909 


5.238 


C in] 


0.299 


2470 


4.049 


[O n] 


0.288 


2697 


3.708 


Ho I 


0.226 


2723 


3.672 


He I 


0.222 


2764 


3.618 


He I 


0.217 


2829 


3.535 


He I 


0.209 


2945 


3.396 


He I 


0.199 


3188 


3.137 


He I 


0.195 


3697 


2.705 


H I 


0.163 


3704 


2.700 


H I 


0.162 


3712 


2.694 


H I 


0.161 


3722 


2.687 


H I 


0.160 


3729 


2.682 


[O II] 


0.160 


3734 


2.678 


H I 


0.159 


3750 


2.667 


H I 


0.158 


3771 


2.652 


H I 


0.155 


3798 


2.633 


H I 


0.152 


3820 


2.618 


He I 


0.150 


3835 


2.608 


H I 


0.148 


3889 


2.571 


He I 


0.142 


3965 


2.522 


He I 


0.134 


3970 


2.519 


H I 


0.133 


4026 


2.484 


He I 


0.126 


4102 


2.438 


H I 


0.117 


4340 


2.304 


H I 


0.086 


4363 


2.292 


[O ml 


0.082 


4388 


2.279 


He I 


0.078 


4400 


2.273 


Bband 


0.076 


4471 


2.237 


He I 


0.064 


4713 


2.122 


He I 


0.024 


4861 


2.057 


H I 


0.000 


4922 


2.032 


He I 


-0.009 


4959 


2.017 


[O in] 


-0.015 


5007 


1.997 


[O in] 


-0.022 


5016 


1.994 


He I 


-0.024 


5500 


1.818 


Vband 


-0.091 


5518 


1.812 


[CI III] 


-0.093 


5538 


1.806 


[CI III] 


-0.096 


5755 


1.738 


[N II] 


-0.123 


5876 


1.702 


He I 


-0.138 


6548 


1.527 


[N II] 


-0.218 


6563 


1.524 


H I 


-0.220 


6583 


1.519 


[N II] 


-0.222 


6678 


1.497 


He I 


-0.233 


6716 


1.489 


[S ii] 


-0.238 


6731 


1.486 


[S ii] 


-0.239 


7065 


1.415 


He I 


-0.278 


7136 


1.401 


[Ar III] 


-0.286 


7236 


1.382 


C ii 


-0.298 


7254 


1.379 


O i 


-0.300 


7281 


1.373 


He 1 


-0.303 


7319 


1.366 


[O II] 


-0.307 


7330 


1.364 


[O II] 


-0.309 


9229 


1.084 


H I 


-0.507 


9464 


1.057 


He I 


-0.529 


9546 


1.048 


H I 


-0.537 


9603 


1.041 


He I 


-0.542 


10049 


0.995 


H I 


-0.579 


10311 


0.970 


He I 


-0.598 
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Table 3. Observed ratios of common-upper-level lines, hni/hmzi and predictions from 

theory 



Dataset 


Theory 


Uncorrected 


Corrected for reddening 




Z82 a 


W96 b 




using H I lines 


FOS-1SW 


0.75 


0.81 


0.35 ±0.01 


0.71 ±0.01 


STIS-SLITlb 


0.75 


0.81 


0.19 ±0.01 


0.61 ±0.02 


STIS-SLITlc 


0.75 


0.81 


0.27 ±0.01 


0.66 ±0.01 



a Using Zeippen (1982) transition probabilities 
b Using Wiese et al. (1996) transition probabilities 
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Table 4. He I relative intensities, Ix/Iun, after correcting for reddening using H I Balmer 

(and Paschen) lines and our extinction curve 



A 


n 




Predicted 








Observed 






1 A ) 




Case B 


Model M 


Mn It' 


OTV92 




BVV00 


EPG04 


PPIQ 1 QW 


o ± lo-oi_il 1 ±c 








(A \ 




{*>) 


K<) 


(8) 


(9) 


(1U) 


fill 












2 


3 S — n 3 P 










3889 


3 


2.315 


0.780 


0.580 


1.019 zb 0.362 


1.165 z_t 0.242 


1.126 ± 0.191 


0.950 ± 0.152 


1.005 ± 0.131 




+H I 




4.575 a 


3.389 


2.651 


3.112 ± 0.561 


3.518 ± 0.393 


3.702 ± 0.163 


3.267 ± 0.103 


3.425 ± 0.080 




3188 


4 


0.878 


0.441 


0.337 


0.933 ± 0.487 






0.536 ± 0.043 


0.521 ± 0.018 




2945 


5 


0.414 


0.290 


0.229 










0.310 ± 0.017 


0.257 ± 0.010 


2829 


6 


0.228 


0.200 


0.162 














+ [Fe iv] 






0.268 


0.248 










0.232 ± 0.017 




2764 


7 


0.141 


0.140 


0.116 










0.140 ± 0.030 


0.120 ± 0.004 


2723 


8 


0.093 


0.100 


0.085 










0.093 ± 0.028 


0.085 ± 0.003 


2696 


9 


0.065 


0.073 


0.063 










0.048 ± 0.026 














2 


l S - n l P 










5016 


3 


0.566 


0.517 


0.488 




0.496 ± 0.070 


0.551 ± 0.014 


0.502 ± 0.007 




0.473 ± 0.013 


3965 


4 


0.221 


0.206 


0. 193 






0.236 ± 0.006 


0.209 ± 0.007 






-t-H l/[Ne in] 






4.973 


4.187 


4.405 i 0.794 








4 994 + n 117 




3614 




0.108 


0.104 


0.095 


0.146 ± 0.037 


0.085 ± 0.019 


0.109 ± 0.003 


0.104 ± 0.007 






3448 


6 


0.059 


0.061 


0.052 


0.105 ± 0.055 






0.068 ± 0.006 






3355 


7 


0.037 


0.041 


0.032 


0.208 ± 0.109 






0.042 ± 0.006 






3297 


g 


0.024 


0.029 


0.021 








0.027 ± 0.008 
















2 


3 P — n J S 










7065 


3 


0.612 


1.714 


1.833 


1.414 ± 0.255 




1.951 ± 0.079 


1.832 ± 0.130 


1.617 ± 0.065 


1.458 ± 0.024 


4713 


4 


0.103 


0.153 


0.157 


0.150 =t 0.038 


0.145 ± 0.032 


0.153 ± 0.005 


0.151 ± 0.002 






4121 




0.038 


0.046 


0.046 


0.046 ± 0.013 


0.049 ± 0.011 


0.044 ± 0.001 


0.049 ± 0.002 






3868 


6 


0.018 


0.020 


0.020 






0.019 ± 0.003 


0.017 ± 0.002 






_j-[Ne m] 






2.846 


2.846 


3.382 ± 0.610 












3733 


7 


0.010 


0.011 


0.011 






0.010 ± 0.001 


0.011 ± 0.004 






3652 


8 


0.006 


0.007 


0.007 








0.005 ± 0.002 






3599 


9 


0.004 


0.004 


0.004 














3563 


10 


0.003 


0.003 


0.003 














3537 


11 


0.002 


0.002 


0.002 








0.003 ± 0.001 
















2 


3 P - n S D 










5876 


3 


2.789 


2.815 


2.820 


2.736 ± 0.493 


3.119 ± 0.349 


3.255 ± 0.121 


3.165 ± 0.100 


2.731 ± 0.058 


2.902 ± 0.044 


4471 


4 


1.000 


1.000 


1.000 


1.000 ± 0.150 


1.000 ± 0.100 


1.000 ± 0.023 


1.000 ± 0.010 


1.000 ± 0.019 


1.000 ± 0.011 


4026 


5 


0.472 


0.466 


0.466 


0.484 ± 0.103 


0.440 ± 0.062 


0.497 ± 0.014 


0.479 ± 0.015 


0.431 ± 0.021 




3820 


6 


0.254 


0.254 


0.254 


0.237 ± 0.059 


0.240 ± 0.054 


0.232 ± 0.006 


0.262 ± 0.008 


0.227 ± 0.027 




3705 


7 


0.155 


0.155 


0.155 


Blend 




0.176 ± 0.005 


0.153 ± 0.008 






3634 


8 


0.102 


0.102 


0.102 


0.129 ± 0.038 


0.096 ± 0.021 


0.129 ± 0.004 


0.104 ± 0.007 






3587 


9 


0.071 


0.071 


0.071 


0.090 ± 0.026 


0.075 ± 0.017 


0.083 ± 0.003 


0.071 ± 0.006 






3554 


10 


0.052 


0.051 


0.051 


0.071 ± 0.021 


0.038 ± 0.008 


0.057 ± 0.002 


0.050 ± 0.005 






3531 


11 


0.038 


0.038 


0.037 


0.060 ± 0.017 




0.037 ± 0.002 


0.039 ± 0.007 






3513 


12 


0.029 


0.029 


0.029 


0.046 ± 0.013 




0.034 ± 0.002 


0.028 ± 0.005 






3499 


13 


0.023 


0.023 


0.023 


0.042 ± 0.022 




0.030 ± 0.002 


0.023 ± 0.005 






3488 


14 


0.019 


0.018 


0.018 


0.028 ± 0.014 






0.018 ± 0.004 






3479 


15 


0.015 


0.015 


0.015 


0.015 ± 0.008 






0.013 ± 0.003 






3474 


16 


0.012 


0.013 


0.012 








0.013 ± 0.004 






3468 


17 


0.010 


0.011 


0.011 








0.007 ± 0.003 
















2 


'P - n L S 










7281 


3 


0.151 


0.155 


0.145 


0.129 ± 0.032 




0.153 ± 0.004 


0.152 ± 0.012 




0.148 ± 0.002 


5048 


4 


0.038 


0.040 


0.036 




0.042 ± 0.009 


0.043 ± 0.002 


0.126 ± 0.003 






4438 


5 


0.015 


0.016 


0.015 


0.017 ± 0.006 


0.013 ± 0.004 


0.0149 ± 0.0004 


0.016 ± 0.001 






4169 


6 


0.008 


0.008 


0.007 


0.009 ± 0.004 


0.011 ± 0.003 


0.0111 ± 0.0008 


0.013 ± 0.001 






4024 


7 


0.005 


0.005 


0.004 






0.0055 ± 0.0005 
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Table 4 — Continued 



A n Predicted Observed 

(A) Case B Model M Model K OTV92 EPTE98 BVVOO EPG04 FOS-1SW STIS-SLITlc 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



2 1 P - n 1 D 



6678 


3 


0.784 


0.737 


0.731 


0.704 ± 0.149 


0.939 ± 0.133 


0.955 ± 0.024 


0.912 ± 0.055 


0.729 ± 0.031 


0.743 ± 0.013 


4922 


4 


0.272 


0.258 


0.256 


0.238 ± 0.051 


0.263 ± 0.037 


0.289 ± 0.007 


0.267 ± 0.004 


0.265 ± 0.028 


0.204 ± 0.005 


4388 


5 


0.125 


0.120 


0.119 


0.112 ± 0.033 


0.117 ± 0.026 


0.121 ± 0.003 


0.120 ± 0.003 


0.103 ± 0.016 


0.105 ± 0.003 


4144 


6 


0.067 


0.066 


0.065 


0.058 ± 0.017 


0.060 ± 0.013 


0.067 ± 0.002 


0.063 ± 0.003 






4009 


7 


0.040 


0.040 


0.040 


0.042 ± 0.012 


0.046 ± 0.010 


0.040 ± 0.001 


0.037 ± 0.002 






3927 


8 


0.026 


0.026 


0.026 


























3 "S - n 










9464 


5 


0.023 


0.053 


0.054 


0.080 ± 0.020 






0.027 ± 0.004 






8362 


6 


0.015 


0.030 


0.029 


Blend 






0.033 ± 0.004 






7816 


7 


0.010 


0.018 


0.017 


0.014 ± 0.004 






0.022 ± 0.002 






7500 


8 


0.007 


0.012 


0.011 


0.012 ± 0.004 






0.014 ± 0.001 






7298 


9 


0.005 


0.008 


0.007 


0.008 ± 0.003 




0.0090 =t 0.0003 


0.0095 ± 0.0010 






7161 


10 


0.004 


0.006 


0.005 






0.0067 ± 0.0002 


0.0070 ± 0.0007 






7062 


11 


0.003 


0.004 


0.004 






0.0046 ± 0.0002 


0.0048 ± 0.0005 






6989 


12 


0.002 


0.003 


0.003 






0.0030 ± 0.0001 


0.0032 ± 0.0004 






6934 


13 


0.002 


0.003 


0.002 






0.0025 ± 0.0002 


0.0033 ± 0.0005 






6890 


14 


0.001 


0.002 


0.002 














6856 


15 


0.001 


0.002 


0.002 








0.0022 ± 0.0004 


















3 1 S - n 










9603 


6 


0.005 


0.006 


0.005 


0.011 ± 0.003 












8915 


7 


0.004 


0.004 


0.003 


0.007 ± 0.002 






0.0057 ± 0.0009 






8518 


8 


0.003 


0.003 


0.002 


0.004 ± 0.002 






0.0029 ± 0.0005 






8266 


9 


0.002 


0.002 


0.002 


Blend 












8094 


10 


0.001 


0.002 


0.001 


0.005 ± 0.002 






0.0015 ± 0.0003 






7972 


11 


0.001 


0.002 










0.0012 ± 0.0003 


















3 S P - n 










10311 


6 


0.029 


0.029 


0.029 


0.012 ± 0.006 






0.038 ± 0.006 






9517 


7 


0.019 


0.019 


0.019 








0.009 ± 0.001 






9063 


8 


0.013 


0.013 


0.013 








0.015 ± 0.002 






8777 


9 


0.009 


0.009 


0.009 


0.010 ± 0.003 






0.024 ± 0.003 






8583 


10 


0.007 


0.007 


0.007 














8445 


11 


0.005 


0.005 


0.005 














8342 


12 


0.004 


0.004 


0.004 








0.0067 ± 0.0009 






8265 


13 


0.003 


0.003 


0.003 


Blend 












8204 


14 


0.003 


0.003 


0.003 








0.0026 ± 0.0004 






8156 


15 


0.002 


0.002 


0.002 








0.0022 ± 0.0004 






8116 


16 


0.002 


0.002 


0.002 








0.0016 ± 0.0003 






8084 


17 


0.001 


0.002 


0.002 








0.0007 ± 0.0003 






C Hf 3 (f 


rom 


H I lines) 






0.61 ± 0.11 


0.50 ± 0.21 


0.30 ± 0.12 


0.82 ± 0.04 


0.52 ± 0.05 


0.66 ± 0.04 


ChB (from 


H I + He I 


subset) 




0.57 ± 0.06 


0.60 ± 0.12 


0.31 ± 0.07 


0.82 ± 0.07 


0.52 ± 0.03 


0.69 ± 0.05 


He + /H 


+ (f 


rom H I + He I subset) 




0.098 ± 0.007 


0.084 ± 0.007 


0.081 ± 0.005 


0.088 ± 0.003 


0.083 ± 0.002 


0.095 ± 0.003 


Hc+/H 


+ (f 


rom source 


select He I 


lines) 


0.089 ± 0.002 


0.089 ± 0.009 




0.087 ± 0.001 







a Case B H I + He I blend determined using a typical He+ /H+ (0.088) as found from these data 
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Table 5. CLOUDY input parameters 



Quantity 


Model M a 


Model K b 


T cS (K) 


35200 


41200 


log0(H) 


13.05 


13.10 


radius (pc) 


0.27 


0.27 


logn H (inner) (cm' 3 ) 


3.4 


3.4 


Turbulence (km s _1 ) 


12 


12 



a Mihalas (1972) stellar atmosphere 
b Kurucz (1979) stellar atmosphere 
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Table 6. Model abundances relative to H (12+log(X/H)) 



Element 


Model M a 


Model K b 


(1) 


(2) 


(3) 


He 


10.98 


10.98 


C 


8.48 


8.42 


N 


7.85 


7.73 





8.60 


8.65 


Ne 


7.78 


8.05 


S 


7.00 


7.22 


Ar 


6.48 


6.62 


CI 


5.00 


5.46 


Fe 


6.48 


6.48 



a CLOUDY H II region abun- 
dances from Baldwin et al. (1991); 
Rubin et al. (1991) and OTV92 

b EPG04 abundances 
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Table 7. Constraints on Model Parameters 



Quantity 


Orion Nebula (1SW) 




SLITlc Model M Model K 


SB(H/3) a 


66.8 ±0.7 


64.1 


63.2 


A5007/H/3 


3.13 ±0.05 


3.5 


4.3 


A5007/A3726 b 


4.0 ±0.1 


5.9 


6.6 


(A4959 + A5007)/A4363 c 


389 ±7 


421.0 


392.1 


(A6548 + A6583)/A5755 c 


64 ± 1 


61.6 


65.0 


A6731/A6716 d 


2.00 ± 0.04 


2.1 


2.0 


A3726/A3729 d 


2.3 ±0.2 


2.5 


2.4 


A1907/A1909 d 


1.20 ±0.04 


1.1 


1.2 


(A3726 + A3729)/A7325 d 


6.4 ±0.6 


4.7 


6.0 



a Surface brightness in units of 10~ 13 ergs cm -2 s _1 arcsec 

ionization indicator 

c Electron temperature indicator 

d Electron density indicator 



-32 - 



Table 8. Dereddened HST STIS V2 and FOS observations (i "a / 'In/3, predicted) 



ID 


A (A) 




SLIT1 


3 




SLIT1 






SLIT 


2 


3 




SLFI 


2 






FOS-1SW 


C in] 


1907 


o 


0691 


± 





0038 


o 


0872 


± 
































C in] 






0591 


± 







j? 


0726 


± 































[U IIJ 


2470° 




0718 


± 





nn'« 




1208 


± 









0550 


± 




0024 




0807 


± 




0032 




lOlO It 


He I 


2697° 












































nn9i -4- n nn i 9 


He I 


2723 c 


o 


0036 


± 





0078 


o 


0042 


± 





0002 






















Q 


0041 ± 0.0012 









0075 


± 





0007 





0059 


± 





0002 
























nn^i 4- n nn i 'i 


ric [x 1 c iv J 


9899 












































m no -t n nnnv 
U±Uz ± u.uuu / 




2945 




0126 


± 









0126 


± 





























m -t n nnnv 
U±oD ± u.uuu / 


He I 


3188 










































Q 


Ui/y it u.uuu* 


H I 


369 7 C 


o 


0161 


± 





0326 


o 


0089 


± 





0006 


{1 


0175 


± 


o 


0140 


o 


0061 


± 


o 


0123 












0152 


± 









0187 


± 







jj 


0256 


± 








0194 


± 










H I 


3712° 




0153 


± 





0017 




0182 


± 





nnnu 




0138 


± 




0171 




0139 


± 




0042 






xj T _i_rq tttI 

STL 1— p 1 O XXX J 


3722 


() 


0298 


± 





0015 


n 


0363 


± 





0011 




0305 


± 




0010 




0370 


± 




0025 






ro ttI 

11J 


3726 


o 


6725 


± 





0199 


o 


7875 


± 





0267 


o 


6690 


± 


o 


0187 


o 


7199 


± 


o 


0330 








3729 


o 


3234 


± 





0298 


o 


3451 


± 





0298 


o 


3353 


± 


o 


0252 


o 


3503 


± 


o 


0199 






H I 


3734 c 


o 


0297 


± 





0027 


o 


0266 


± 





0013 


o 


0201 


± 


o 


0007 


o 


0247 


± 


o 


0018 






H I 


3750 c 


o 


0265 


± 





0022 


o 


0306 


± 





0011 


o 


0330 


± 


() 


0028 


o 


0260 


± 


() 


0013 


o 


0285 it 0.0018 


II I 


3771 c 


o 


0319 


± 





0015 


o 


0387 


± 





0011 


o 


0380 


± 


() 


0014 


o 


0377 


± 


() 


0015 


o 


0397 it 0.0014 


H I 


3798° 





0425 


± 





0011 





0509 


± 





0014 





0400 


± 





0012 





0490 


± 





0019 





0520 it 0.0015 


Ho I 


3820° 













































0099 it 0.0012 


He I+H I 


3889 













































1499 ± 0.0027 


He I+H l+[Ne III] 


3965 













































2155 ± 0.0041 


He I 


4026° 













































0188 ± 0.0009 


H I 


4340° 





4817 


± 





0058 





4817 


± 





0126 





4542 


± 





0042 





4822 


± 





0158 





4784 ± 0.0065 


[O III] 


4363 





0109 


± 





0003 





0109 


± 





0003 





0088 


± 





0003 





0136 


± 





0005 






He I 


4388° 





0057 


± 





0002 





0051 


± 





0002 





0046 


± 





0002 





0057 


± 





0002 





0045 ± 0.0007 


He I 


4471 c 





0457 


± 





0006 





0481 


± 





0013 





0415 


± 





0005 





0452 


± 





0015 





0436 ± 0.0010 


H I 


4861 c 


1 


0085 


± 





0112 


1 


0044 


± 





0260 


1 


0382 


± 





0076 





9990 


± 





0323 





9671 ± 0.0146 


He I 


4922° 





0092 


± 





0010 





0097 


± 





0003 





0145 


± 





0010 





0104 


± 





0007 





0115 ± 0.0012 


[O III] 


4959 


1 


2109 


± 





0145 


1 


0535 


± 





0272 


1 


2355 


± 





0097 


1 


1675 


± 





0380 






[O III] 


5007 


3 


6142 


± 





0393 


3 


1408 


± 





0809 


3 


6795 


± 





0269 


3 


4824 


± 





1124 






He I 


5016 





0260 


± 





0007 





0226 


± 





0008 





0213 


± 





0006 





0196 


± 





0007 






[CI in] 


5518 





0040 


± 





0002 





0033 


± 





0001 





0025 


± 





0002 





0030 


± 





0002 






[CI in] 


5538 





0062 


± 





0002 





0055 


± 





0001 





0043 


± 





0001 





0052 


± 





0002 






[N II] 


5755 





0074 


± 





0002 





0122 


± 





0003 





0068 


± 





0002 





0084 


± 





0003 






He I 


5876° 





1363 


± 





0015 





1365 


± 





0035 





1345 


± 





0008 





1371 


± 





0044 





1179 ± 0.0019 


[N II] 


6548 





1507 


± 





0020 





1877 


± 





0049 





1501 


± 





0015 





1587 


± 





0052 






H I 


6563° 


2 


8274 


± 





0327 


2 


7863 


± 





0741 


2 


8477 


± 





0209 


2 


8694 


± 





0949 


3 


0933 it 0.0466 


H I 


6563° 










































2 


7514 ± 0.0402 


[N II] 


6583 





4601 


± 





0065 





5834 


± 





0153 





4582 


± 





0050 





4896 


± 





0160 






He I 


6678° 





0343 


± 





0005 





0346 


± 





0009 





0332 


± 





0003 





0337 


± 





0012 





0317 ± 0.0023 


He I 


6678 c 













































0313 ± 0.0012 


[S ii] 


6716 





0213 


± 





0005 





0217 


± 





0006 





0247 


± 





0005 





0209 


± 





0007 






[S ii] 


6731 





0392 


± 





0007 





0432 


± 





0012 





0455 


± 





0006 





0409 


± 





0014 






He I 


7065 





0669 


± 





0008 





0675 


± 





0018 





0668 


± 





0005 





0668 


± 





0022 





0694 ± 0.0026 


[Ar III] 


7136 





1499 


± 





0019 




1547 


± 





0041 





1460 


± 





0013 




1542 


± 





0051 






C ii 


7236 





0026 


± 





0002 




0026 


± 





0002 





0024 


± 





0001 




0023 


± 





0002 






O i 


7254 





0012 


± 





0001 




0011 


± 





0000 





0009 


± 





0001 




0012 


± 





0001 






He I 


7281 c 





0065 


± 





0001 




0068 


± 





0002 





0058 


± 





0001 




0061 


± 





0002 






[O II] 







0613 







0008 




0942 







0025 





0549 







0006 




0655 







0022 





2203 ± 0.0047 


[O ii] 


7330° 





0500 


± 





0006 




0775 


± 





0020 





0442 


± 





0005 




0541 


± 





0018 




Blend 



C H 3 0.89 ± 0.05 0.71 ±0.04 0.74 ± 0.08 0.87 ± 0.06 0.54 ± 0.02 

He+/H+ 0.090 ± 0.003 0.094 ± 0.003 0.089 ± 0.004 0.089 ± 0.003 0.083 ± 0.002 



a FOS/G780 grating 
b FOS/G570 grating 

c incs used to determine reddening curve. 
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Fig. 1.— HST STIS Slit 1 (52" x 0.5", right) and Slit 2 (left) shown overlapping with HST 
FOS positions 1SW (right circle) and x2 (left circle) (Baldwin et al. 1996; Rubin et al. 1997), 
respectively. Slit positions are overlaid on F656N (Ha) WFPC2 image. 9 1 Ori C is marked 
with ax. N is up and E to the left. 
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Ha H/S H6 




2 3 4 5 6 7 

wavenumber C/im -1 ) 



Fig. 2. — Analytic stellar-based extinction curves, (top) Dash-dotted line: CCM89 extinction 
curve presented here in differential form, E\_y / Eb-v {Rv = 5-5) for ease of comparison with 
original BS81 data (x). BS81 data should be adjusted by roughly +0.25 mag - the result 
of a calibration error pointed out by CC88, whose corrected stellar data are also plotted 
(*). Solid line: our modification of the CCM89 formulation to better fit the CC88 corrected 
data and shape near 2.7 //m _1 . {bottom) Same as above, but the window has been scaled to 
accentuate the differences between CCM89 and our modified curve. 
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Fig. 3. — CP70 extinction curve (dashed line) including extension into the UV using un- 
corrected BS81 data following Walter (1993). Our new extinction curve (solid curve) is 
overplotted, showing an unexplained discrepancy in the IR portion of the CP70 curve. The 
renormalized CP70 nebular extinction curve (dotted line) shows a marked improvement to 
the IR discrepancy. The UV portion of the scaled CP70 curve is rederived directly from the 
BS81 Ex^v/Eb-v data, using the new scaled jy and fs- 
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8000 6000 
wavelength (A) 



4000 




10' 



8000 6000 
wavelength (A) 



4000 



Fig. 4. — Best-fit differential extinction curves (dashed line) as determined from EPG04 H I 
Balmer and Paschen lines. Residuals and error bars are plotted for each H I emission line (+) 
used in the fit. (top) our modified CCM89 extinction curve, Ry = 5.5, Cup = 0.82 ± 0.04. 
Also shown are data from common upper-level pairs of [S n] (•) and He I triplets (A) 
and singlets (A) after applying this reddening correction. The optical member of the pair, 
usually with the lower error, has been plotted as I cor rected = {predicted- (bottom) Residuals 
for extinction correction according to our scaled CP70 nebular curve, Cup = 0.76 ± 0.04. 
The common-upper level pairs are not included, but would result in a similar scatter about 

log(-^ corrected/ ' I predicted) 0. 
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Fig. 5. — Energy-level diagram of He I, showing (top) singlet and (bottom) triplet transitions. 
For the sake of clarity, not all observed transitions are shown. 
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5 10 15 20 25 

T (X3889) 

Fig. 6. — Variation in line flux (I\/Iuji) as a function of To(A3889), the optical depth at the 
center of A3889 for v tur b = (solid line), (top) /7065//4471; (bottom) Isssg/Iun- For the same 
physical models, the changes when v tU rb = 5, 10, 20, and 30 km s -1 are also shown. 
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Fig. 7. — Variation of He I lines (2 3 S' — n 3 P series) as a function of v turb in the full 
nebular model. The hatched boxes represent the dereddened FOS-1SW observations (and 
uncertainties) of AA7065, 3889 (deblended), 3188, and 2945 (from Table 4). 
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Fig. 8.— Observed He I decrements (OTV92, •; BVVOO, *; EPG04, ■; FOS-1SW, A; STIS- 
SLITlc, ♦) compared to case B (open symbols) and CLOUDY model M (filled symbols) 
predictions for four series. The corrected values, I C orrected, are determined from our extinction 
curve, which is calibrated for each dataset using H I and a subset of He I lines (see § 4.5). The 
68.3% confidence intervals are shown for each member, n, of the series. For visual clarity, 
each set of observed data has been offset by 10A, and the uncertainties have been excluded 



for the OTV92 data in the 2 3 P 



n 



series (in the UV, these uncertainties are close to 



±0.2 dex). Note the change in scale between the top and bottom figures. 
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7000 6000 5000 4000 3000 
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Fig. 9. — Differential extinction curve for FOS-1SW (dashed line) presented as C\ — Ch/3- 
Residuals, log(I correc t e d/ ^predicted) (relative to case B), and la errors are plotted for each 
emission line used in the fit: H I (+), He I (triplets, A; singlets, A) and [On] (x). The three 
He I lines most affected by radiative transfer (AA7065, 3188, and 2945) are not included in 
the fit, but are plotted relative to case B values (o) and model M predictions (•). 
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Fig. 10. — Like Fig. 9, but for STIS-SLITlc and a different scale. The two He I lines 
most affected by radiative transfer (AA7065 and 2945) are not included in the fit, but are 
plotted relative to case B values (o) and model M predictions (•). For case B, A7065 has 
\og(I correc ted I Ipredicted)= 0.4. Similarly shown is singlet A5016 which was also not included in 
the fit. 



